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ABSTRACT 

We focus on a modified version of Horava - Lifschitz theory and, in particular, we 
consider tlie impact of its weak- field static spherically symmetric limit on the galaxy 
dynamics. In a previous paper, we used the modified gravitational potential obtained 
in this theory to evaluate the Milky Way rotation curve using a spheroidal truncated 
power - law bulge and a double exponential disc as the only sources of the gravitational 
field and showed that the modified rotation curved is not in agreement with the data. 
Making a step forward, we here include also the contribution from a dark matter halo 
in order to sec whether this helps fitting the rotation curve data. As a test case, we 
consider a sample of spiral galaxies with smooth baryon matter distribution and well 
measured circular velocity profiles. It turns out that, although a marginal agreement 
with the data can be found, this can only be obtained if the dark matter halo has an 
unrealistically small virial mass and incredibly large concentration. Such results can 
be interpreted as a strong evidence against the reliability of the gravitational potential 
obtained in the modified version of Horava- Lifschitz theory that we consider. 

Key words: gravitation - dark matter ~ galaxies : kinematic and dynamics 



1 INTRODUCTION 

In an attempt to construct a theory of quantum gravity 
in 3 -I- 1 dimensions, decoupled from string theory, Ho- 
rava proposed, in 2009, a theory which took advantage of 
theoretical conc epts investigat ed in the field of condensed 
matter physics (|Horaval l2009al ') . Precisely, Horava's theory 
is analogue to a scalar -field model developed some years 
ago by Lifschitz in the framework of the theory of quan- 
tum critical phenomena, in which the full Lorentz symme- 
try emerges only at infrared (IR) fixed point. For this rea- 
son, one usually refers to it as the Horava - Lifschitz (HL) 
theory. The main feature of this theory as a field quan- 
tum gravity model is the anisotropic scaling of the space 
and time dimensions. Because of this anisotropic scaling, 
one finds the time dimension plays a privileged role show- 
ing up in a reduced invariance: full difFeomorphism invari- 
ance is abandoned, while only a subset, foliation - preserving 
diffeomorphism, is kept. This foliation has the property of 
becoming coincident with the standard diffeomorphism in- 
variance i n the IR limit, where General Relativity (GR) is 
recovered (|Horavall2009bl ) . By the way, the properties caus- 



ing the enthusiastic reception it has gained among physi- 
cists are certainly its be ing power - cou nting renormaliz- 
able and definitely finite (|Horavall2009d ). Many of i ts fea- 
tures , like its ultraviolet (UV) a nd IR behaviour (jVisserl 
[200^, its impact on cosmo logy jCalcagni 1 12009|) . spher- 



ically symmetric solutions (ILUj Mei fc Popel |2009| ). black 
holes (Kehagias fc Sfctsos 2009) an d So l ar Sy stem tests 
(|Harko et al.ll201ll ; ilorio fc Rug gicra,2oTol . l201lh have been 
deeply analyzed. Unfortunately, its original formulation, 
based on two main assumptions, that is the detailed bal- 
ance condition (that reduces the number of operators in the 
action and simplifies some properties of the quantum sys- 
tem) and the projectabiUty condition (a certain part of the 
spacetime metric, the so-called lapse function, can be set 
globally t o unity), soon s howed lots of problems both in the 
IR limits (| Orlando fc Reff crt 2009) and intrinsic pari ty vio- 
lation in the purely gravitational sector of the mode (|Visseil 
I2OO9I ). The attempt to fix those problems led to abandon the 
detailed balan ce condition, still retaining the projectabil- 
ity c ondition (^Kehagias & Sfctsos, '2009|; iKiristis fc Kofina^ 
I2OIOI : iCapasso fc Polvchronakos 1,20101 ). Although the 



2 V.F. Cardone et al. 



discussion about the foundations and the possible phe- 
nomenological and conceptual implications of the H L theory 
and of its various modifications is still open (see e.g. lSotirioul 
l|201[J ) and references therein), we believe it is worth sys- 
tematically investigating its consequences at every scales. 
Follo wing this path, in a previous paper (| Cardone et al. I 
I2OIOI ) we have complemented recent works on the cosmolog- 
ical consequences of the modified HL gravity as presented 
in Sotiriou et al. (2009a, b) by addressing its impact on the 
gravitational potential : in contrast to the results obtained 
in Tang & Chen (2009), we showed that static spherically 
symmetric solutions other than the Schwartzschild - de Sitter 
one do exist. Precisely, we worked out a general formalism 
to estimate the rotation curve (that is, the circular velocity, 
Vc, as a function of the distance from the centre, R) for an 
extended source, showing that the HL theory can boost the 
Vc{R) with respect to the Newtonian value, as a consequence 
of the modification induced on the gravitational potential of 
a point mass. As an application, we evaluated the Milky 
Way rotation curve using a spheroidal truncated power-law 
bulge and a double exponential disc as the only sources of the 
gravitational field, finding that the modified rotation curve 
is unable to fit the data. This means that the modified HL 
model taken into account can not play the role of an alter- 
native solution to the missing mass problem. Motivated by 
this conclusion, we therefore take here a step further adding 
a dark matter halo to the galaxy modelling. Provided the 
galaxy model is reliable, one should then be able to fit the 
rotation curve of spiral galaxies within the framework of the 
modified HL potential. Moreover, such a test would also al- 
low us to constrain the modified potential parameters and 
hence put some constraints on the coupling quantities enter- 
ing the HL Lagrangian. On the contrary, should the model 
provide a poor fit to the data, one can draw some inter- 
esting lessons on the validity of the HL gravity from a test 
performed on galactic scales thus complementing the theo- 
retical investigations and the observational analyses based 
on cosmological data. It is interesting to point out that even 
though our work is based on a static spherically symmetric 
solution obtained in the framework of HL gravity, our results 
can be applied as well to an arbitrary perturbation of the 
Newtonian gravitational potential having the same form. 

The plan of the paper is as follows. In Sect. 2, we review 
the derivation of the rotation curve for extended systems in 
the framework of Horava- Lifschitz gravity and then apply 
it to a spiral galaxy modelled as the sum of a thin disc and 
spherical dark matter halo. The spiral galaxies sample and 
its main properties are described in Sect. 3 together with 
the method used to constrain both the halo and theory pa- 
rameters. The results are presented and discussed in Sect. 4, 
while summary and conclusions are given in Sect. 5. 



2 HORAVA - LIFSCHITZ ROTATION CURVE 

As stated above, we consider here the version of the HL 
gravity proposed by Sotiriou, Visser & Weinfurtner (2009 a, 
b) whose weak field limit has been fir st discussed by Tang fc 
Chen (2010). In our previous paper (|Cardone et al. ||2010| ). 
we extended the Tang & Chen results by looking for a so- 
lution different than their Schwartzschild - de Sitter one. We 



indeed found that the gravitational potential generated by 
a point mass m read : 

^{r) =^N{r) + ^HL{r) , (1) 
where $jv(r-) — —Gm/r is the Newtonian potential and 



j4_ 2 _B _ D 1 
<i>HL[r) - ^2/'' + 4/ 12/r-i 



(2) 



is the correction due to the HL theory. Here, the (^4, B, D) 
parameters are related to the coupling constants entering the 
HL Lagrangian, while / is an unknown function of the source 
mass m satisfying the condition f(m = 0) = 1. The correc- 
tive t erm may be conveniently rewritten as (| Cardone et al. I 
l2010t ) : 



'i>HL{r) = 



n_\ _ Br^ J]_ 



.(3) 



with Mq the Sun mass and an arbitrary chosen reference 
radius introduced to define the dimensionless quantity 77 = 
r/r^. In Eq.©, we have finally defined the scaling radii : 
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where quantities labelled with a are evaluated for /i = 
m/Mo = 1. It is worth noting that the second term in Eq.Q 
simply adds a constant to the potential which has no effect 
in any situation of interest, so that we will henceforth neglect 
this term. We stress that this by no way means that B can 
be set to zero. Indeed, B — means / = 1 so that also 
D vanishes (again, see I Cardone et al. I (|2010l ) for details) 
and we go back to the Schwartzschild - de Sitter solution, 
while we are here interested in the more general case. We 
therefore assume B ^ 0, but nevertheless hereafter neglect 
its contribution because it drops off from the derivation of 
the quantities we are interested in. 

As Eq.((3l) refers to the case of a point source mass, we 
have first to generalize it to an extended system in order to 
derive the rotation curve Vc{R) = R{d^/dR). To this end, 
we have also to take care of the fact that the corrective term 
$£f_L does depend on m in a way that we do not explicitly 
know, so that a nonlinear dependence can not be excluded 
a priori. As such, the superposition principle does not hold 
anymore and a different procedure has to be worked out . 
This problem has been addressed in ICardone etHI (|2010h . 
where we derived a general expression which holds for any 
potential leading to a point mass gravitational force factor- 
izable as : 



with and fr{r) dimensionless functions depending on the 
particular form of the point mass gravitational potential $p. 
Remembering that Fp — — V$p, it is only a matter of alge- 
bra to show that : 

/m(m) = M , friv) = 1/??^ . 

for the Newtonian potential, while they are : 
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for the firs t and second te r m in ^hl- Following the steps 
detailed in ICardone et all (|20ld ). one finally gets for the 
rotation curve the expression : 



vl (R) = G po Rd r] 



X / r] dr] 
10 



Piv'X')dC' / MAo)d9', (5) 
J Jo 



with 



/21 1/2 



(6) 



Ao = A(6I = C = 0) = [r,^ + J?'^ - 27777' cose' + C'^] 
and i^d-t) mass function (hereafter MF) normalized so that 

fitl){fi)dp = Po/Mq , 



with Po the mass density in the neighborhood of the refer- 
ence radius. 

It is worth stressing that Eq.© is fully general and can 
be used to compute the rotation curve, provided the expres- 
sion for ffiip) and /r(77) are given. As a consistency check, it 
is easy to show that, for the Newtonian potential, the term 
depending on the MF is identically equal to unity, so that 
Eq.© reduces to a simple rewriting of the standard result. 
For the HL term, we get a dependence on the MF through 
the multiplicative term on the second row of Eq.(0. Note 
that the MF here only plays the role of scaling up or down 
the rotation curve. This is a consequence of the HL correc- 
tive term not depending on m. As such, one has simply to 
sum the contribution of all the stars notwithstanding their 
mass and this is indeed what the MF term gives in the HL 
case. The total rotation curve for the HL theory will finally 
be given as : 

vl[R) = vlriR) + vl^iR) = vl[R) + v\[R) -f vl{R) , 

where one can use the results known in literature for the 
computation of the Newtonian term v% (R) and our general 
rule with f^{p) and frirf) given above to estimate the HL 
contribution vjj]^{R). Note that, in order to evaluate these 
latter terms, we need to know not only the mass density (as 
in the Newtonian only case), but also the local MF. Since we 
do not know this function, we can not separately constrain 
the parameters (r^o, 7-00). We can, however, easily skip this 
problem as follows. Let us consider the va term and write 
it explicitly as : 



Iaq) 
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We can now define a new length scale 77A as : 

'^^^'^n i::: ^(^'^^ — 1 

so that the first HL corrective term simply becomes : 



(8) 



■"a{R,i1a) = 



G po Rd 77 (2/77" 
1^ dr] 



p{v'X')dC' / Aode' (9) 
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which can be evaluated, as function of 77A, even if the local 
MF is not known. Using a similar trick for the v^{R), we 
can now fit the model to the rotation curve data thus con- 
straining the parameters (77^,7713). Should the local MF be 
known and a theoretically motivated expression for f{p) be 
available, one can then translate the estimates of (rjA,r]D) 
into constraints on (77^10, »7d0). 

Eq.® makes it clear that the shape of the local MF and 
the functional expression of f{p) affects the overall shape of 
the total rotation curve through a rescaling of the corrective 
terms scalelengths (77^,770). On the other hand, it is possi- 
ble that one has to set these latter quantities on a case - by - 
case basis in order to fit the rotation curves data. Actually, 
this is not a problem for the consistency of the theory since 
(77a,?7d) should not be universal quantities being related to 
the HL Lagrangian couplings only through integrals involv- 
ing both f{p) and tp{p). However, since we expect f{p) to 
be the same for all systems, one has to assume that the local 
MF is different from one system to another in order to make 
a non universality of {rjARdjVDRd) consistent with the con- 
stancy of {ta, td)- Actually, whether the MF is universal or 
not is still an open question which we will not address here. 



3 SPIRAL GALAXIES ROTATION CURVES 

Historically, the flatness of spiral galaxies rotation curves 
was the first and (for a long time) more convincing evide nce 
for the existence of dark matter (|Sofue &: Rubin ||2001^ . It 
is therefore a mandatory test for a whatever modified grav- 
ity theory to first demonstrate it is able to reproduce the 
observed circular velocity profile. Some care is, however, 
needed to be sure that the test results are not biased by 
possible systematics in both the theoretical prediction and 
the observational measurements. First, it is worth stressing 
that Eq.(l5]) implicitly assumes that the main contribution 
to Vc{R) comes from stars moving on circular orbits in the 
disk plane. On the contrary, the asymmetric drift and the 
contribution from stars moving in the spiral arms may lead 
to non circular motions which are not taken into account 
in the above derivation. Moreover, since the HL corrective 
terms have to be computed numerically, a clumpy matter 
distribution can introduce an artificial bias related to how 
the density profile has been sampled and interpolated. In or- 
der to avoid both these problems, we therefore have to select 
spiral galaxies dominated by the stellar component and with 
a small contribution from the spiral arms such as Sbc, Sc, 
Sd systems, checking that the disk surface brightness profile 
can indeed be well approximated by the exponential model 
we have considered above. 

From an observational point of view, the desiderata for 
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the rotation curve data are mainly dictated by the need to 
sample the circular velocity profile with both high precision 
(so that significant constraints can be derived) and over a 
large radial range (in order to probe deep into the outer halo 
dominated regions). It is also highly preferable that the disk 
inclination has been well determined so that no systematic 
bicis is induced on the measured Vc{R) values. Finally, it is 
worth remembering that what one actually measures is the 
circular velocity as a function of the radial distance from the 
galaxy centre measured in arcsec. In order to convert to the 
linear units to get Vc — Vc{r) with r in kpc, we need to know 
wit great accuracy the galaxy distance from us. An uncer- 
tainty on this quantity translates in a global rescaling of the 
mass related parameters thus making difficult to raw any 
definitive conclusion on the importance of the dark matter 
contribution to the total mass budget. 

Fulfilling all the above constraints is not an easy task 
and greatly reduce the sample of candidates for our analy- 
sis. Actually, it is worth stressing that our aim here is not 
to constrain the theory parameters, but rather to check its 
overall viability. From this point of view, a single Sc galaxy 
with an exponential disk and rotation curve data well sam- 
pled out to large radius is preferable to a large sample of 
systems with poorly determined rotation curves and domi- 
nated by not smooth disks. Indeed, should the HL modified 
potential fail to reproduce the data for this golden probe, 
one could safely argue against the theory as a whole. 

3.1 The data 

Motivated by the above considerat ions, we have looked at 
the THINGS (|de Blok et al. 1 120081 ) dataset comprising 19 
galaxies with high quality rotation curves data finely sam- 
pling a large radial range (up to ~ lORd). Both the inclina- 
tion and the distance of these systems have been determined 
with high precision, while a careful analysis of the non cir- 
cular motions have shown that these latter are negligible for 
most of the cases. The THINGS galaxies are, therefore, ideal 
candidates for our analysis fulfilling all the observational re- 
quirements we have hinted at in the above discussion. 

Notwithstanding their high quality circular velocity 
measurements, not all the galaxies in the THINGS sample 
can be used for our analysis. Indeed, in order to not in- 
troduce spurious errors in the the theoretical rotation curve 
derivation, we ask that the stellar components have a smooth 
profile. To this regard, it is worth remembering that the lu- 
minosity profiles of the THINGS galaxies are typically de- 
scribed as the sum of an exponential disk and a second inner 
disk referred to as the bulge (only present for some galax- 
ies). The mass density profile is then obtained multiplying 
the surface brightness by a stellar M/L ratio inferred from 
the galaxy color. Since the color can radially change, the 
M/L ratio could also not be constant so that the final den- 
sity can be significantly different from the classical exponen- 
tial profile and present undesirable wiggles and ripples. In 
order to avoid possible problems with the interpolation of 
non smooth components, we therefore select only galaxies 
with a smooth component and set the total disk mass Md 
and scalelength radius Rd in such a way that the Newtonian 
rotation curve reported in the THINGS data is the same as 
the one of a Freeman (1970) disk (see later) with the same 
{Md,Rd) values. Such a smoothness cut gives a sample of 



Table 1. Properties of sample galaxies. Explanation of the 
columns : name of the galaxy; logarithm of the bulge total mass 
(in solar units); bulge scalelength (in kpc); logarithm of the disk 
total mass (in solar units); disk scalelength (in kpc). Values for 
the bulge are not reported if there is no evidence for such a com- 
ponent in the surface brightness profile. 



Id 


log Mb 




log Md 


Rd 


NGC 2841 


10.41 


0.67 


11.04 


4.13 


NGC 3621 






10.24 


2.20 


NGC 5055 


9.33 


0.53 


11.08 


3.71 


NGC 6946 


9.55 


0.45 


10.78 


3.24 



only six galaxies. After excluding two systems with evidences 
for strong irregularities in the circular velocity profile, we se- 
lect the four galaxies in Table[T] briefly presented below. 

• NGC 2841. This is an early -type (Sb) spiral with 
absolute B mag Afs = —21.21 and distan ce D — 14. 1 Mpc . 
Its rotation curve was first measured by iBeeeman I (jl987[ ) 
and shows a good overall agreement with the THINGS 
data we use here. Similarly, the position angle and the 
inclination differ by only few degrees from P. A . = 153° and 
i = 74° recommended bv lde Blok et al. I (120081 ). 

• NGC 3621. At a distance D = 6.6 Mpc, NGC 3621 is a 
late -type spiral with a regular HI distribution. Although 
it is quite luminous with Mb = —20.06, no other deter- 
mination of the rotation curve was available before the 
THINGS one which have also determine the position angle 
and inclination as 345° and 65°, respectively. 

• NGC 5055. This Sbc galaxy (D = 10.1 Mpc, 
Mb = —21.12) has an extended and warped tenuous 
outer HI disc which originates the wiggle visible for 
large R values in its rotation curve. The THINGS rota- 
tion curve agree s well with the previous determination 
(|Battaglia et al. 1 |2006| ) and extends further out thus 
being able to probe the warp for the first time. Position 
angle {P.A. — 102°) and inclination {i — 59°) are also in 
agreement. 

• NGC 6946. Although well known for its population of 
HI structures, the dynamics of NGC 6946, a Ms = -20.61 
late -type spiral at D = 5.9 Mpc, has been relatively 
little studied because of its low inclination (i = 33°, 
P.A. = 243°). The offset between th e THINGS data we use 
here and previo us determinations (jCarignan et al. 1 Il990l : 
iBoomsma 1120071 ) is likely due to different assumptions for i. 

Although quite limited in number, such a small sample is 
nevertheless ideal for our test of the HL potential. Indeed, 
the regularity of the mass components, the detailed sam- 
pling of the rotation curve, the unambiguous determination 
of both the inclination and the distance leave almost no 
space for questioning the validity of the fit results. Indeed, 
should the fit work bad, one can confidently argue that it 
is a failure of the theory and not an outcome of some unac- 
counted for theoretical or observational bias. Motivated by 
this consideration, we prefer to not relax our severe selection 
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criteria thus losing statistics, but improving the reliability of 
any inference on the validity of the HL modified potential. 



3.2 Modelling spiral galcixies 

In lCardone et al. I (j201(]| ) , we have attempted to fit the Milky 
Way rotation curve using the modified HL potential and the 
luminous disk, only in order to see whether it is possible to 
give off the dark matter contribution still obtaining a flat 
circular velocity profile. Such a test unambiguously demon- 
strated that this is not the case, so that a dark halo is still 
unavoidable to fit the rotation curve data. Actually, such 
a result could be anticipated looking at the scaling with r 
of the additive terms in the point mass gravitational po- 
tential. Indeed, while the one oc r~* quickly decreases, the 
correction scaling as should be weighted by a very large 
va in order to avoid being dominant in the inner regions 
where the rotation curve is consistent with the prediction 
from Newtonian gravity. On the contrary, ta must not be 
too large if one wants to boost the circular velocity in the 
usually dark matter dominated region. Balancing these two 
opposite desiderata is actually not possible, so that a dark 
halo is still needed to fit the data. We will therefore model 
spiral galaxies as two component systemj^, namely the stars 
(in the bulge and the disk) and the dark matter. 

We assume the stars (both for the bulge and the disk) 
are distributed in an infinitely thin and circularly symmetric 
disk. The surface density profile simply reads : 

AT 

S(i?) = ^exp(-i?/i?,) , 

with {Md,Rd) the total mass and the sc alelength radiu s. 
The Newtonian rotation curve is given by ([Freeman lllQTOl ) : 



«dJv(-R) = 



27rGEoiid(W2) 
[Io{v/2)Ko{v/2) 



■hivmKiiv/2)] , (10) 



with T] = R/Rd and Ki Bessel functions of order I of the 
first and second type, respectively. 

While the observed photometry motivates the use of the 
exponential profile for the disk, the choice of the dark halo 
model is not trivial. Numerical simulations of structure for- 
mation are typically invoked as a direct evidence favouring 



the use of the NFW density law (iNavarro et al. _ _ . _ _ . . 
or its variants (|Moore et al. Il99aljing fc Suto l(2000l '). How- 



:e tavourmg 
19961 . 119971 ) 



ever, the NFW model is the outcome of DM only for sim- 
ulations performed in a Newtonian framework, while here 
we are working in a modified gravity theory. In principle, 
one should therefore rely on the results of simulations which 
include both the effect of the different potential and the im- 
pact on the evolution of structures due to deviations from 
GR. While this has been done in some cases, for in stance in 
the context of f{R) theories (|Schmidt et al. II2OO9I '). we are 
not aware of any similar test for HL gravity. We can never- 
theless make some qualitative arguments to drive the choice 
of the halo profile. Firstly, we note that, in the intermediate 
regions, it is reasonable to expect that the r~* corrective 



^ Note that we are here neglecting the gas component since, for 
the galaxies considered, its contribution to the mass budget is 
much smaller than both the stars and the dark matter. 



term has yet faded away, while the one is still too small 
if rA is large enough. In this same limit, we can also extrap- 
olate that, close to the virial radius, the potential is still 
approximately Newtonian, so that the results from N - body 
simulations can be taken as a first order approximation. On 
the contrary, in the inner regions, the term can give a sig- 
nificant contribution, so that one can not rely on the scaling 
with r inferred from Newtonian simulations, but a different 
density profile could be expected. An efficient way to param- 
eterize these qualitative c onsiderations is rep resented by the 
generalized NFW model jjing fc Suto Il200(t l reading : 



Ph{r) 



with 



4-KR^g{R^ir/Rs 
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r 

R's 



-(3-c) 



2), 



(11) 



(12) 



and B{x,a,b) the Beta function. In Ea. (|lip . Mvir and Rvir 
are the virial mass and radius, respectively. They are not 
independent, as the following relation holds 



Rvi. 



( ZM^r 



1/3 



with Ath the overdensity for spherical collapse and pM = 
SH^ flM/SnG the mean ma tter density today. We fol- 
low [Br^n^^orman] (|l998l ) for Ath and se t (QM ,h) = 
(0.28, 0.70) in accordance with lKomatsu et al. I (|201GI 1. Note 
that the value of Ath has been derived for a concordance 
ACDM model, hence implicitly assuming that GR rather 
than HL is the theory of gravity. However, from the cos- 
mological point of view, HL reduces to ACDM for a spa- 
tially fiat homogeneous and isotropic universe, so that no 
bias is induced by this assumption. Because of the spherical 
symmetry, the Newtonian rotation curve may then be easily 
evaluated as : 



VhN{r) 



GMhjr) ^ GM^ir gjr/Rs) 



(13) 



while the HL corrective term will be computed numerically. 
Note that, to this end, we need to assign the three halo pa- 
rameters which we choose to be the inner regions logarithmic 
density slope, a, the circular velocity at the virial radius, 
= GMvir/Rvir, and the concentration Cyir — Rvir/Rs- 
Note that we prefer to use {Vvir, Cvir) instead of {Mvir,Rs) 
as fitting parameters since these are better constrained by 
the data being related, respectively, to the overall scale of 
the rotation curve and the balance between the stellar and 
dark matter terms in the intermediate radial regions. 



3.3 Fitting procedure 

In order to constrain both the HL and halo parameters, 
we employ a standard Bayesian approach, first defining the 
likelihood function as : 



£(p) (X exp 



exp ■ 



X^(P) 



-IE 



(14) 



6 V.F. Cardone et al. 



Table 2. Best fit solutions and marginalized constraints on the model parameters. Columns are as follows: 1. galaxy id, 2. best fit 
parameters, 3., 4., 5., 6., 7., 8. median and 68% confidence ranges for the fitted quantities. 



Id 






P6/ 






a 


c, 


jir 




ir 


log 


VA 


log riD 


NGC 2841 


(1.82, 


0.21, 


731, 


88, 2.95, 


-0.65) 


1 93+0.21 
^•=''-'-0.19 




396 


+ 328 
-223 


5lj 


-32 
-21 


i.6i: 


fO.84 
-0.47 


—0 66+0-35 


NGC 3621 


(1.05, 


0.08, 


211, 


25, 1.67, 


-1.11) 


1 14+0.23 


17+0-28 


291 


+330 
-190 


25j 


-14 
-14 


0.63; 


+0.39 
-0.17 


n 84+0.25 
-0.84_p 44 


NGC 5055 


(0.71, 


0.08, 


746, 


75, 2.32, 


-0.38) 


0.79^H^ 




168 


+388 
-80 


25+27 
^°-13 


1.14; 


f 1.05 
-0.27 


-0 60+°-^'' 


NGC 6946 


(0.78 


, 0.24 


, 16, 


7, 0.96, - 


1.10) 




17+0-28 


267 


+323 
-203 


26j 


-31 
-IC 


0.72; 


+0.59 
-0.15 


n 77+0.31 
"■' '-0.61 



where p is the set of the parameters of the model, v°'"'{Ri) 
and v^J^{Ri) are the observed (with a measurement error £;) 
and theoretically predicted values of the circular velocity at 
the radius Ri of the i-th point, and the sum runs over the 
observed data points. The best fit is obtained by maximizing 
the likelihood -C(p), but it is worth stressing that, according 
to the Bayesian philosophy, the best estimate of the param- 
eter Pi is not the best fit one. On the contrary, one has to 
marginalize over the remaining parameters and look at the 
shape of the marginalized likelihood function defined as : 

Ci{pi) on j C{p)dpi . . . dpi-idpi+i . . . dpn , 

with n the total number of parameters. Actually, what we 
did was running a Monte Carlo Markov Chain (MCMC) 
code to efficiently explore the five dimensional parameter 
space and use of the histogram of the values for the param- 
eter Pi to estimate the mean, the median and the 68% and 
95% confidence ranges. Note that, because of degeneracies 
among the parameters of the model, the best fit parameters 
Pbf may also differ from the maximum likelihood ones Pml, 
i.e. the set obtained by maximizing each of the marginalized 
likelihood functions. 

It is mandatory to explain which are the parameters 
we constrain by our fitting procedure. Firstly, we scale the 
HL characteristic radii (rA,r_D) with respect to the disc 
scalelength Rd, and then skip to logarithmic units using 
(log rjA , log rjn ) as model parameters in order to explore 
a wider range. As already quoted above, the halo model 
is assigned by three parameters which we choose to be 
(q, Cvir, Vvir), while the virial mass is estimated a posteriori 
from the V^ir value. Finally, although the total disc mass 
has been set from the beginning, it is worth noting that 
such an estimate has been obtained converting colors into 
a fiducial M/L ratio T fid assuming a diet Salpeter initial 
mass function (IMF) and a given recipe for the different in- 
gredients entering the stellar population synthesis code. As 
a consequence, one can not exclude a priori that the actual 
M/L ratio T does not equal the fiducial one so that we add 
/'^ f id to the list of fitting parameters 



4 RESULTS 

In order to constrain the parameters of the model, we have 
run, for each galaxy, three Markov Chains using the Cell- 



man - Rubin (1992) criterium to check the convergence of the 
algorithm. The best fit solution and the constraints on the 
parameters are summarized in Table 2, while Fig.[T] shows 
the best fit curve superimposed to the data. 

As a preliminary remark, we remember the reader that 
the errors Si on each data point are not Gaussian distributed, 
since they also take into account systematic misalignments 
between HI and Ha measurements and other effects lead- 
ing to a conservative overestimate of the true uncertainties. 
The final error on each point is, indeed, computed by adding 
in quadrature the dispersion around the best fit tilted ring 
model estimate and a pseudo Icr uncertainty due to the dif- 
ferences between the appro aching and receding s ides of the 
circular velocity profile (see Ide Blok et al. |[2008l for further 
details). Note that such procedure leads to larger errors for 
galaxies with residual non circular motions and/or warped 
HI distribution in the outer regions. As a consequence of 
the non Gaussian uncertainties, a value = /dof ~ 1 
for the best fit model (with dof = TV — 6 the number of 
degrees of freedom) could also be an evidence for an overes- 
timate of the uncertainties rather than the outcome of the 
model closely matching the data. 

That this is indeed the case may be easily understood 
noting that we get = 0.92 for NGC 2841 which is clearly 
not well fitted by the model. On the contrary, for the other 
systems, is a reasonably good indicator being — 1.05 
for NGC 3621, = 3.62 for NGC 5055, x^ = 1.74 for NGC 
6469, in qualitative agreement with what one expects from 
the overall agreement between the model and the data. In 
order to get some insight on the agreement between the 
model prediction and the data, it is worth briefiy discussing 
each galaxy on a case -by -case basis. 

NGC 2841- This is the worst fitted galaxy in the sam- 
ple with only the very inner regions well fitted by the data. 
While the data shows a decreasing trend, the best fit circu- 
lar velocity profile starts diverging after an initial decrease. 
Notwithstanding the relatively small value, we can nev- 
ertheless safely argue that the modified potential is not able 
to fit the NGC 2841 rotation curve. 

As a further evidence, we can also consider the 
marginalized constraints on the model parameters. First, we 
note that T /T fid takes quite large values. Since the diet 
Salpeter IMF used to infer T fid provides the highest possi- 
ble M/L value compatible with the TuUy- Fisher relation, it 
is quite hard to conceive a particular combination of stellar 
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Figure 1. Best fit curves superimposed to the data for NGC 2841 (top left), NGC 3621 (top right), NGC 5055 (bottom left), NGC 6469 
(bottom right). See Table 2 for the values of the best fit parameters. 



population properties able to give values of T/T fid ^ 1.25 
so that we can safely deem as unrealistic the T/T fid values 
needed to reconcile the model with the data. 

Similarly, the concentration parameter takes extremely 
large values. It is worth noting, however, that what we are 
referring to as concentration is not the usually defined one. 
Indeed, one typically refers to = Rvir/R-2 with i?_2 
the radius where the logarithmic slope of the density pro- 
files equals -2. For the gNFW model we are considering, it 
is C-2 = c^ir/(2 — a) so that, for the median {oi,Cvir) val- 
ues in Table 2, we get C-2 ~ 200 which is still quite large 
if compared to the values in literature (c_2 ~ 10 — 15). 
Although previous estimates have been obtained in a New- 
tonian framework, it is unlikely that such strongly concen- 
trated haloes are indeed realistic. It is therefore worth won- 
dering what is driving the fit towards such large c^ir values. 
To this end, let us first note that the estimated virial mass 
IS given bylogM,„ = 10.79t°j* which is typical for a spiral 
galaxy. Since c^tr — Rvir/Rs and Rvir takes typical values, 
we must conclude that it is Ra to be unrealistically small. 

In order to understand why such values are preferred, 
we note that the smaller is Rs, the earlier the regime 
r/Rs » 1 is achieved and the halo density profile starts 
scaling as . As a consequence, for R » Rb, the 
Newtonian contribution to the halo tends to quickly vanish 
and the outer regions circular velocity is mainly determined 
by the HL term. On the contrary, in the very inner 
regions, R/Rs « 1, the halo Newtonian Vc{R) becomes 



quite small and its role is played by the HL r"** term. 
Indeed, the values of (log 77A, log t^d) are tailored in such 
a way to compensate a missing halo term in the outer 
and inner regions, respectively. Should Rs not have been 
extremely small, the halo Newtonian term would contribute 
significantly in these two regions thus leading to a mismatch 
with the data. 

NGC 3621. Contrary to NGC 2841, this galaxy is weU 
fitted by the modified potential. However, a deeper look at 
the data shows that, while the inner regions are perfectly 
matched, the agreement between the theoretical curve and 
the data, although overall good, is getting worse as R in- 
creases. Indeed, while the data suggests a fiat behaviour, 
the best fit circular velocity starts increasing linearly as a 
consequence of the term in the HL potential leading to 
Vc{R) (X R for R/Rd » rjA- Should one be able to probe 
Vc{R) to larger R, the model will likely fail to fit the ex- 
tended dataset leading to a situation similar to NGC 2841. 

Table 2 shows that the concentration is still quite large 
so that the same qualitative discussion made for NGC 
2841 also holds here. Moreover, the estimated virial mass, 
log M„ir = 9.7815 03) turns out to be smaller than the disc 
mass which can lead to serious problems on cosmological 
scales (see later for a discussion of this issue). 

NGC 5055. The striking feature of this fit is the oppo- 
site behaviour of the theoretical and observational rotation 
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curve in the very inner regions. Although such a strong dis- 
agreement for a best fit model could be surprising, it is worth 
noting that actually only two points are missed so that their 
weight in the overall is quite small compared to the other 
ones probing the intermediate and outer regions. What we 
have suggested for NGC 3621 clearly takes place for this 
gala^xy. The rotation curve data stays roughly flat (or slowly 
decreasing at most), while the theoretical profile first de- 
creases and then starts increasing linearly as consequence of 
the set in of the term in the IfL modified potential. 

Again, we find a very large concentration indicating 
an extremely small Ra so that the qualitative discussion 
made for NGC 2841 holds for this case too. Moreover, the 
virial mass, logAf^ir = 9.84tQ g^, is one order of magnitude 
smaller than the disc mass so that, contrary to the Newto- 
nian case, the system is stars dominated everywhere. 

NGC 6469. From an observational point of view, NGC 
6469 can be considered a classical example of a galaxy with 
a rotation curve first increasing almost linearly and then 
setting in a long flat profile (although with some wiggles 
likely due to a clumpy gas distribution). On the contrary, 
the best fit curve poorly matches the rising part and shows 
a definitely increasing trend which does not follow the trend 
probed by the data. Moreover, we again find large concentra- 
tions and a virial mass, logMi„> = 9.92liJ'25, smaller than 
the disc mass. 

4.1 Discussion 

What the above results tell us is that our disks -I- halo model 
gives rise to a theoretical rotation curve which can not al- 
ways be reconciled with the observed one. While the theo- 
retical best fit circular velocity provide a satisfactory agree- 
ment with the data for NGC 3621 and NGC 5055 (exclud- 
ing the very inner points), it traces the data for NGC 2841 
and NGC 6469 quite poorly. Moreover, in all cases, there is 
marked trend indicating a rising rotation curve in the outer 
regions, whereas an extrapolation of the data suggests a fiat 
profile. Although one could consider this latter argument 
questionable, it is worth stressing that, in order the fit to 
be acceptable, the price to pay is to have large concentra- 
tions and small halo masses. Being the rotation curves used 
of high quality and not affected by any systematic due to 
wrong inclination or distance determination, we can safely 
argue that the poor matches among the best fit theoretical 
curve and the measured one is a strong evidence of a failure 
of the modified HL potential. 

Two further considerations strengthen this conclusion. 
First, we have qualitatively explained the large concentra- 
tions as a consequence of the need to reduce as much as pos- 
sible the halo contribution in the very inner and very outer 
regions. For the median values of logr^A and log rjo, the two 
corrective terms introduced by the HL modified potential 
boost the Newtonian circular velocity for R/Rd » rjA and 
R/Rd « rjD- Since the disc parameters are held fixed, the 
halo scaling quantities (mass and radius) have to be adjusted 
in order to not lead to Vc{R) values larger than the observed 
ones. In particular, this drives the best fit model towards 
halo virial masses smaller than the disc ones. Considering 
the large error bars, one can at most have log M^ir — log Md, 
which, if extrapolated to the cosmological scales, would im- 



ply Q,b — ^CDM, quite at odds with what is expected from 
the expansion and growth of structures dat41. 

Table 2 shows that the 68% confidence ranges of 
(log r^A, log 77d) for the different galaxies overlap quite well 
so that, although a larger sample is needed, one could ar- 
gue that the scaled radii (77a,»7d) are universal quantities. 
Actually, such a result would represent a shortcoming for 
the consistency of the theory. Indeed, since Rd is different 
from one galaxy to another, one should postulate a coupling 
between the properties of the galaxy and the coupling pa- 
rameters entering the HL Lagrangian (and determining the 
(f A, td) radii) leading to a universal value of (r]A, rjo)- This 
would lead to a gravity theory depending on the particular 
system one is considering which is clearly at odds with the 
fundamental properties the gravity Lagrangian must have. 
A possible way out would be to invoke a non universality of 
the local MF, but it is hard to find a physically motivated 
recipe linking the disc scale radius Rd and the MF in such 
a finely tuned way to make {r]A,r]D) universal quantities. 

Motivated by the above analysis, we can therefore com- 
plement the result in Paper I (where no halo was considered) 
concluding that the HL modified gravitational potential we 
have focused on is unable to reproduce the spiral galaxies 
rotation curve because of a too large boost to the circular 
velocity of both the visible and dark matter components. 



5 CONCLUSIONS 

Initially motivated by its attractive features from the point 
of view of quantum gravity, the HL proposal has soon be- 
come one of the most investigated theories of gravity alter- 
native to GR. In I Cardone et alTI (|2010l ). we have comple- 
mented recent works on its cosmological consequences by 
addressing its impact on the gravitational potential finding 
that, under some suitable conditions on the couplings con- 
stants entering its Lagrangian, it is possible to work out a 
modified gravitational potential made out of the Newtonian 
1/r term corrected by the addition of two further terms scal- 
ing as and r~*, respectively. We have here made a step 
further investigating the viability of this modified potential 
from the point of view of galaxy dynamics. To this end, we 
have tried fitting the rotation curve data of 4 spiral galax- 
ies including both visible (bulge and disc) and dark matter 
components. As a result, we have found that the theoretical 
rotation curve may match the observed one only in the very 
inner regions, but the theory parameters {rjA,rjD) can not 
be adjusted to fit also the outer regions. Moreover, the halo 
virial mass and concentration should take unreasonably odd 
values, the halo mass and characteristic radius being smaller 
than their disc counterparts. While one could look at such a 
low dark matter content as a good news, it is worth stressing 
that HL gravity still needs substantial dark matter on cos- 
mological scales in order to provide the correct evolution for 
the scale factor. As such, lowering the galaxies dark matter 
content does not solve the missing mass problem, but sim- 
ply turns it into the problem of understanding where dark 

^ It is worth stressing that, for a spatially flat universe, the Fried- 
mann equations for the HL gravity reduce to those for the concor- 
dance ACDM model so that one expects the same matter content 
in both theories in order to flt, e.g., the SNela Hubble diagram. 
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matter should be present in order to fill the gap with the 
needed cosmological amount. 

It is worth stressing that, although the above results 
have been obtained assuming a particular halo model, we 
expect that they are qualitatively model independent. As 
discussed above, lowering the halo mass and scale radius is 
needed in order to compensate the boost in circular velocity 
due to the and terms of the HL modified poten- 
tial. Since they dominate in the inner and outer regions, 
respectively, one has to reduce the total mass in these re- 
gions to not override the observed Vc{R). Having fixed the 
stars content, the only way to achieve to do this is indeed to 
have a large concentration and a small virial mass. Chang- 
ing the halo model will likely change the actual values of 
{cjjir, Mvir), but can not reduce the boost due to the HL 
corrective terms. As a consequence, we can safely argue that 
similar results should finally be obtained thus qualitatively 
confirming the analysis presented here. 

As a final remark, we point out that our discussion 
suggests that also the Schwartzschild - de Sitter solution 
(which presents the term, but not the r~'^ one), obtained 
in the framework of SVW version of HL gravity, can 
hardly fit the galaxies rotation curves. We can consequently 
conclude that the analysis of galaxies dynamics represents 
a strong evidence against the reliability of any gravitational 
potential obtained in this modified version of HL gravity. 
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